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摘要 ”蜜蜂 的 级 型 分 化 被 证 实 是 由 蜂王 浆 中 的 Royalactin 决 定 ， 工 蜂 和 蜂王 幼虫 在 级 型 分 化 时 编码 基因 也 
研究 ， 我 们 发 现 ， 在 蜜蜂 幼虫 的 级 型 分 化 过 程 中 ，IncRNA 也 有 着 显著 的 表达 差异 ， 因 此 认为 ，IncRN 的 级 型 
分 化 过 程 . 进一步 的 分 析 显 示 ，lncRNA 可 能 通过 影响 上 下 游 基 因 的 转录 和 功能 执行 的 方式 ， 在 蜜蜂 早期 发 育 的 多 细胞 组 织 
发 育 、 神 经 系统 发 育 和 转录 调控 的 过 程 中 起 到 重要 的 调控 作用 . 
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关键 词 ”蜜蜂 ， 级 型 分 化 ，IncRNA 
学 科 分 类 号 Q75, Q96, S89 DOI: 10.16476/j.pibb.2015.0110 


雌性 蜜蜂 的 级 型 分 化 是 一 种 典型 的 发 育 可 塑性 。” ”转录 本 表达 情况 ， 对 其 中 IncRNA 的 表达 进行 了 分 
现象 点， 同样 的 受精 卵 在 不 同 的 环境 下 接受 不 同 的 析 ， 从 而 找 出 IncRNA 在 蜜蜂 级 型 分 化 过 程 中 的 
食物 导致 了 完全 不 同 的 两 种 发 育 途 径 ， 蜂王 具有 更 ”作用 . 
大 的 体型 ， 更 长 的 寿 俞 ， 完 整 的 生殖 能 力 ， 在 蜂巢 ad Ens 
三 中 负责 交配 和 产 卵 ， 工 蜂 体型 较 小 ， 寿 命 短 ， 不 其 L E 
P 有 生殖 能 力 ， 在 蜂巢 中 执行 各 种 职能 ， 包 括 筑 造 蜂 — L1 单 雄 受精 蜂 群 构建 


f. BW RR. MEIER PRA. REJE 实验 采用 “ 浙 农 大 1 号 ”意大利 蜜蜂 (Apis 
在 发 育 过 程 中 一 直 以 蜂王 浆 为 食物 ， 工 蜂 幼虫 只 在 ”mellifera)， 由 浙江 大 学 负责 孵育 和 收集 ， 为 了 尽 可 


发 育 早期 得 到 了 少量 的 工蜂 浆 作为 食物 ， 之 后 都 是 ”能 减 小 遗传 背景 的 影响 ， 采 用 单 雄 受精 的 蜂 群 进行 
以 蜂蜜 为 食物 .蜂王 浆 中 富 含 多 种 和 蛋白质 、 糖 类 和 实验， 这 种 蜂 群 中 的 蜂王 通过 人 工 受 精 ， 只 接受 一 
核酸 P77， 其 中 的 一 种 和 蛋白质 一 一 Royalactin 被 认为 。 只 雄 蜂 的 精子 ， 从 而 保证 整个 蜂 群 中 所 有 后 代 均 为 
在 蜜蜂 的 级 型 分 化 中 起 到 了 决定 性 的 作用 甸 ， 将 。 同 父 同 母 . 

Royalactin 加 入 到 高 温 灭 活 的 蜂王 浆 中 ， 蜂 王浆 可 12 幼虫 培育 

以 重新 获得 使 幼虫 发 育 为 蜂王 的 能 力 ， 此 外 ， 蜂 了 为 了 收集 同步 化 的 峰 王 和 工蜂 幼虫 ， 让 蜂王 在 
浆 中 所 包含 的 miRNA 也 被 证 实在 蜜蜂 的 级 型 分 化 衬 的 框架 中 24h 内 完成 产 卵 ， 然 后 ， 把 这 个 框架 
过 程 中 发 挥 了 一 定 的 作用 已 ， 对 蜜蜂 幼虫 在 级 型 分 。 转移 到 一 个 空 的 干净 蜂箱 中 .66h 后 ， 转 移 150 只 
此 过 程 中 各 种 基因 的 表达 也 有 过 较 多 的 研究 ， 一 些 
重要 的 基因 在 工蜂 和 蜂王 幼虫 中 有 非常 明显 的 表达 * 中 国 科 学 院 战 略 性 先导 科技 专项 (XDA01020402) 和 国家 高 技术 研 
差异 ， 并 影响 了 工蜂 和 蜂王 不 同 的 发 育 方向 uca. iig (201244020402, 2012AA02A202) 资 助 项 目 . 
研究 报道 ， 长 链 非 编码 RNA(long noncoding +» 通讯 联系 人 

RNA，lncRNA) 在 发 育 和 调控 中 发 挥 了 重要 的 作 E-mail: slchen@zju.edu.cn 

9， 因 此 ， 我 们 对 4—6 日 龄 的 工蜂 和 蜂王 幼 pge Tei oinedesses, bat cnnsibpacen 

虫 进行 了 高 通 量 测 序 ， 检 测 了 所 有 带 poly A 端的 。” 收 稿 日 期 ， 2015-04-22， 接 受 日 期 ， 2015-06-01 
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蜜蜂 幼虫 到 人 造 王 台 ， 同 样 将 转移 后 的 王 台 放 回 原 
来 的 蜂箱 ， 在 第 4 CUI LO a 73~ 90 hb), 第 5 
天 (幼虫 孵化 后 97~ 114), 6 KS) E IRAE Js 
121 ~ 138 分 别 将 普通 框架 和 王 台 中 的 工蜂 和 蜂 
王 幼虫 取出 ， 放 入 液 氮 中 保存 直到 使 用 . 
1.3 总 RNA 提取 ， 文 库 构建 和 高 通 量 测序 

总 RNA 用 Trizol 方法 提取 . W5 pg 总 RNA 
于 构建 文库 . mRNA 文库 按照 Juminas 
mRNA-Seq 8 sample kit 说 明 构 建 ， 先 将 带 poly-A 
端的 RNA 分 离 出 来 ， 然 后 打 断 成 短片 段 ， 反 转 成 
DNA 加 接头 后 通过 PCR 扩 增 构建 成 文库 . 构建 好 
的 文库 通过 Ilumina Hiseq2000 平台 完成 测序 . 
1.4 RNA-seq 数据 分 析 
测序 数据 通过 cufflinks 流程 计算 FPKM 来 检 
测 转录 本 的 表达 水 平 ， 首先 通过 Tophat2 对 测序 所 
得 的 原始 序列 与 蜜蜂 的 基因 组 文件 进行 比 对 ， 能 够 
回帖 到 基因 组 的 序列 用 cufflink 进行 拼接 ， 将 每 个 
样品 中 拼接 好 的 转录 组 整合 成 一 份 基因 组 注释 文 
件 ， 以 此 注释 文件 为 基础 ， 计 算 每 个 转录 本 的 
FPKM (fragments per kilobase per million)， 得 到 的 
结果 通过 cuffnorm 进行 归 一 化 处 理 ， 所 得 结果 作 
为 最 终 的 基因 表达 值 ， 用 来 进行 统计 分 析 . 
1.5 [lJ IncRNA 分 析 

将 差异 表达 的 IncRNA 与 Noncode 数据 库 中 人 


SS 


(a) 


和 小 鼠 的 IncRNA 序列 进行 Blast 比 对 ， 和 筛选 出 其 
中 e-valuec1x10-, bit score>60 的 匹配 作为 同 源 的 
IncRNA， 根 据 人 鼠 中 的 IncRNA 功能 预测 结果 ， 
对 工蜂 幼虫 和 蜂王 幼虫 中 有 表达 差异 的 IncRNA 进 
行 了 功能 预测 . 

1.6 上 下 游 编 码 基 因 的 分 析 

] closetbed 检测 差异 表达 的 IncRNA 上 下 游 
10k 的 编码 基因 ， 用 “reciprocal best hit” 方 法 找 
到 了 蜜蜂 在 果 蝇 中 的 同 源 基因 ， 对 这 些 差 民 表达 的 
IncRNA 上 下 游 基 因 在 果 晶 中 的 同 源 基 因 ， 通 过 普 
林 斯 顿 大 学 Lewis-Sigler Institute 的 Gene Ontology 
Tools 进行 GO 分 析 . 


2 ”结果 与 分 析 


2.1 lncRNA 的 整体 表达 情况 

最 新 的 蜜蜂 基因 组 (Amel 4.5) 中 注释 为 
IncRNA 的 共有 3 993 个 mcm， 在 我 们 的 实验 中 ， 能 
够 检测 到 的 有 565 个 ， 其 中 ， 蜂 王 幼虫 在 4 日 龄 、 
5 日 龄 、6 日 龄 时 可 以 检测 到 的 IncRNA 分 别 有 
427. 378. 336 个 ， 工 蜂 幼 虫 可 以 检测 到 的 
IncRNA 分 别 有 385、362、374 个 . 分别 对 4 HER 
5 日 龄 、6 日 龄 时 工蜂 和 蜂王 幼虫 中 能 够 检测 到 的 
IncRNA 的 表达 进行 对 比 (图 1). 
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Fig. 1 Expression of IncRNA in the queen larvae and worker larvae 


(a) Comparison of IncRNA expression between queen and worker larvae on the 4th day, W4 for worker larvae, Q4 for queen larvae. (b) Comparison of 


IncRNA expression between queen and worker larvae on the 5th day, W5 for worker larvae, Q5 for queen larvae. (c) Comparison of IncRNA expression 


between queen and worker larvae on the 6th day, W6 for worker larvae, Q6 for queen larvae. 


1TE4 Hit. 5 Hit. 6 日 龄 时 都 能 检测 到 的 
IncRNA 有 154 A(X 1). 763% 154 个 共同 表达 的 
IncRNA 中 ，4 日 龄 时 ， 蜂 王 与 工蜂 幼虫 中 表达 差 
异 达 到 2 倍 以 上 的 IncRNA di 66 ^^, 5 日 龄 时 ， 


蜂王 与 工蜂 幼虫 中 表达 差异 达到 2 倍 以 上 的 
IncRNA 有 84 个 ，6 日 龄 时 ， 蜂 王 与 工蜂 幼虫 中 表 
达 差 异 达 到 2 倍 以 上 的 InceRNA 有 114 个 . 在 3 天 
中 的 表达 差异 都 达到 2 倍 的 有 28 个 (图 2). 
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Table 1 The number of identified IncRNAs for all the samples 


Q4 Q5 Q6 WA WS W6 Share 
IncRNA 427 378 336 385 362 374 154 
Coding genes 11274 10 629 10 450 11034 10615 10249 9 307 


*Q for queen larvae, W for worker larvae, the number for age of the larvae. 
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Fig. 2 Distribution of differentially expressed IncRNAs 


我 们 对 所 有 能 够 检测 到 表达 的 IncRNA 的 表达 
量 进行 了 聚 类 分 析 ， 结 果 显 示 ，4 HARIT, TER 
蜂王 的 差异 最 小 ， 随 着 时 间 的 变化 ， 差 异 越 来 越 
K, 6 日 龄 时 ， 表 达 差 异 达 到 最 大 ， 其 中 6 日 龄 的 
工蜂 与 其 他 幼虫 的 差异 最 明显 (图 3). 
2.2 IncRNA 保守 性 分 析 

把 蜜蜂 中 的 IncRNA 和 noncode 数据 库 吧 中 人 
和 小 鼠 中 已 知 的 IncRNA 进行 比 对 ， 找 到 了 蜜蜂 中 
保守 性 相对 比较 高 的 546 个 hcRNA. 这 些 相 对 保 
守 的 IncRNA 在 4 日 龄 、5 日 龄 、6 日 龄 时 表达 差 
异 达 到 2 倍 以 上 的 分 别 有 12、13、28 个 ， 对 这 些 Fig.3 Cluster of IncRNAs for all samples 
IncRNA 在 人 和 小 鼠 中 可 能 的 同 源 IncRNA 的 功能 
预测 显示 ，4 日 龄 时 ，lncRNA 的 功能 主要 分 布 在 
惫 相 和 蛋白 有 反应、 生长 调控 、 记 忆 等 方面 (图 4a). 记忆 等 方面 (图 4b). 6 HARI, ncRNA 的 功能 主 
5 日 龄 时 ， IncRNA 的 功能 主要 分 布 在 信号 转 导 、 要 分 布 在 RNA RAM I ESE. DNA 为 模板 的 转录 
DNA 为 模板 的 转录 调控 、G 蛋白 偶 联 受 体 通路 和 ” ”和 调控 以 及 信号 转 导 等 方面 (图 4c). 
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Fig. 4 Functional prediction of differentially expressed ncRNAs 
The X axis is the number of IncRNA predicted to having the function of this GO term. 


23 ”上 下 游 编码 基因 分 析 

我 们 对 有 表达 差异 的 IncRNA 上 下 游 10k 的 编 
码 基 因 进 行 了 检测 . 4 日 龄 、5 日 龄 、56 日 龄 时 表 
达 差 异 达到 2 倍 的 基因 里 ， 分 别 有 46、55、79 个 
可 以 在 上 下 游 10k 内 找到 编码 基因 . 

用 “reciprocal best hit" 7; 25:79 ?»3 48] T 2X ETE 
果 量 中 所 有 的 同 源 基因 .对 这 些 在 差异 IncRNA 上 


育 、 


调控 主要 集中 在 神经 系统 发 育 、 多 细胞 组 织 发 
转录 调控 等 方面 (图 Sc). 
2.4 编码 基因 的 表达 情况 

除了 IncRNA， 还 检测 了 编码 基因 的 表达 
水 平 ， 最 新 的 蜜蜂 基因 组 注释 OGSv3.2 中 共有 
15314 个 编码 基因 四， 我 们 的 实验 中 共 检 测 到 了 
10566 个 其中， 蜂王 幼虫 在 4 日 龄 、5 日 龄 、6 


下 游 10k 范围 内 找到 的 基因 ， 找 出 其 在 果 晶 中 的 同 
源 基因 并 进行 了 GO ahr. GO 的 结果 显示 ， 这 些 
基因 大 部 分 都 是 调控 基因 . 在 4 日 龄 时 ， 调 控 主 要 
集中 在 多 细胞 组 织 发 育 、 细 胞 分 化 、 多 细胞 组 织 过 
程 、 神 经 系统 发 育 等 方面 (图 Sa). 在 S 日 龄 时 ， 调 
控 主 要 集中 在 信号 、RNA 聚合 酶 开 转 录 、 细 胞 交 
流 、 多 细胞 组 织 过 程 等 方面 (图 5b)， 在 6 日 龄 时 ， 


日 龄 时 可 以 检测 到 的 编码 基因 分 别 有 10 212、 
9855. 9761 个 ， 工 蜂 幼 虫 在 4 日 龄 、5 日 龄 、6 
日 龄 时 可 以 检测 到 的 编码 基因 分 别 有 10 092、 
9819, 9555 个 ， 分 别 对 4 日 龄 、5 日 龄 、6 日 龄 
时 工蜂 和 蜂王 幼虫 中 能 够 检测 到 的 编码 基因 的 表达 
进行 了 对 比 (图 6). 
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Fig. 5 GO analysis of nearby coding genes of differentially expressed IncRNAs 
The X axis is -Ig(P-value). 
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Fig. 6 Expression of coding genes in the queen larvae and worker larvae 
(a) Comparison of coding genes expression between queen and worker larvae on the 4th day, W4 for worker larvae, Q4 for queen larvae. (b) 
Comparison of coding genes expression between queen and worker larvae on the 5th day, W5 for worker larvae, Q5 for queen larvae. (c) Comparison of 


coding genes expression between queen and worker larvae on the 6th day, W6 for worker larvae, Q6 for queen larvae. 
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在 4 日 龄 、5 日 龄 、6 日 龄 时 都 能 检测 到 的 纺 l l 
码 基因 有 8986 个 ， 在 这 些 共同 表达 的 编码 基因 3 vl 论 


中 ，4 日 龄 时 ， 蜂 王 与 工蜂 幼虫 中 表达 差异 达到 2 在 近来 的 研究 中 ， 越 来 越 多 地 发 现 LncRNA 
倍 以 上 的 编码 基因 有 1 083 个 ， 其 中 工蜂 幼虫 中 表 在 名 种 生命 活动 中 发 挥 着 重要 功能 ， 包 括 生长 发 


达 高 的 有 448 个 ， 蜂 王 幼 虫 中 表达 较 高 的 有 635 育 中 习 、 表 观 遗 传 、 癌 症 的 发 生发 展区 、 干 细胞 


个 ，5 日 龄 时 ， 蜂 王 与 工蜂 幼虫 中 表达 差异 达到 2 > 在 高 等 
们 以 上 的 编码 基因 有 2 871 个 ， 其 中 工蜂 幼虫 中 表 动物 中 所 做 的 研究 ， 但 是 越 来 越 多 的 实验 证 明 ， 不 
达 高 的 有 876 TREE SRDS RASTA 1 995 — 只 在 疹 椎 动物 中 ， 在 无 消 椎 动物 中 ，IncRNA 同样 
^ 6 日 龄 时 ， 蜂 王 与 工蜂 幼虫 中 表达 差异 达到 2 也 发 挥 了 重要 的 生物 学 功能 

倍 以 上 的 编码 基因 有 6 744 个 ， 其 中 工蜂 幼虫 中 表 密 峰 是 常用 的 社会 学 和 行为 学 模式 生物 ， 密 蜂 
达 高 的 有 1335 个 ， 蜂 王 幼 虫 中 表达 较 高 的 有 5 409 让 的 级 型 分 化 现象 也 被 广泛 地 关注 和 研究 ， 蜂 王浆 
个 (图 7)， 可 以 看 到 ， 随 着 分 化 的 时 间 越 入 ， 工 蜂 中 的 Royalactin 被 证 实在 蜜蜂 的 级 型 分 化 中 起 到 了 
和 蜂王 幼虫 中 有 表达 差异 的 编码 基因 数目 越 来 越 。 六 定性 作用 ， 工 峰 和 蜂王 幼虫 发 育 中 的 重要 基因 
Z. 但 是 ， 蜂 王 幼 虫 中 表达 量 较 高 的 编码 基因 数目 和 通路 也 有 较 多 看 


始终 多 于 工蜂 幼虫 ， 这 与 之 前 的 研究 结果 -一致 m TE euet. 
miRNA 在 蜜蜂 的 级 型 分 化 中 也 起 到 了 辅助 调控 的 

作用 中 
6000} 在 本 实验 中 ， 通 过 mRNA-seq 方法 ， 检 测 了 工 
5000 - 蜂 和 蜂王 幼虫 中 所 有 带 poly-A 端 的 转录 本 的 表达 
4000| 量 . 虽然 IncRNA 通常 都 带 有 poly-A， 但 是 在 我 们 
3000 L 的 实验 中 ， 仅 有 565 个 可 以 检测 到 有 表达 ， 占 已 注 
PM FE IncRNA 的 比例 只 有 14% ， 而 能 够 检测 到 的 
Al mRNA 占 已 注释 mRNA 的 比例 可 以 达到 66%. 我 
| m 图 | 们 另外 的 实验 中 发 现 ， 蜜 蜂 脑 组 织 中 有 1 196 个 
iu di sa IncRNA 可 以 检测 到 有 表达 .对 这 两 部 分 取 并 集 发 
现 ， 有 335 个 IncRNA 在 两 份 数 据 中 都 可 以 检测 
Fig. 7 The one of genes having higher expression 到 ， 推 测 这 可 能 是 由 于 IncRNA 组 织 特异 和 时 间 特 

只 的 直达 引起 的 

差异 表达 的 IncRNA 随 着 幼虫 的 发 育 而 越 来 越 


多 ， 聚 类 的 结果 与 此 一 致 ， 对 这 些 差 异 表达 的 

IncRNA 进行 了 进一步 分 析 . Noncode 数据 库 吧 中 

我 们 检测 了 IncRNA 上 下 游 10k 以 内 编码 基因 ”收集 了 人 和 小 鼠 中 已 知 的 IncRNA 并 且 对 这 些 

的 表达 情况 ，4 日 龄 、5 Hát. 6 日 龄 时 ， 有 表达 ”IncRNA 进行 了 功能 预测 ， 我 们 把 蜜蜂 中 差异 表达 
差异 的 编码 基因 自身 在 工蜂 和 蜂王 幼虫 中 表达 差异 ”的 IncRNA 与 Noncode 数据 库 中 的 IncRNA 进行 了 

达到 2 倍 的 分 别 有 11、23、45 个 ， 比 例 分 别 达 到 比 对 ， 找 到 其 中 相对 保守 的 IncRNA 及 其 预测 的 功 
T 23.91%、41.82% 和 56.96%( 表 2). 能 .结果 表明 ， 在 4 日 龄 、5 Hie. 6 日 龄 时 ， 

IncRNA 始终 在 信号 传导 、 记 忆 、 生 长 调控 、 转 录 


调控 等 方面 发 挥 作 用 ,在 5 日 龄 、6 日 龄 时 ， 还 在 
Table 2 The expression of nearby genes 多 细胞 组 织 发 育 、DNA 修复 、 有 蛋白 质 修 饰 等 方面 
of differentially expressed IncRNAs RIETEN. 

—-—— 很 多 时 候 ，IncRNA 都 是 通过 调控 上 下 游 的 编 

Diff iall d IncRNA. 64 84 114 $ ae s 
Do EEG 码 基因 来 发 挥 作用 的 ， 所 以 我 们 对 工蜂 和 蜂王 幼虫 
Nearby mRNAs 46 55 79 中 差 异 表达 的 IncRNA ER 10k 的 因 进 和 VN 
Differentially expressed mRNAs 11 23 45 d me ui mi H2 


Percent 23.91% 41.8296 56.96% pr. 结果 显示 ， 绝 大 多 数 基因 都 是 调控 相关 基因 . 
调控 主要 发 生 在 多 细胞 组 织 发 育 、 神 经 系统 发 育 、 
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RNA 转录 酶 开户 动 的 转录 等 方面 ， 而 这 些 基 因 中 


达 


n 


UA 


IncRNA íh 


执 
发 


[8 


[9 


很 多 本 身 在 工蜂 和 蜂王 幼虫 中 就 有 显 阔 的 表达 差 
， 因 此 推测 ，IncRNA 很 可 能 通过 影响 这 些 上 下 
基因 的 转录 和 功能 执行 来 发 挥 调控 作用 . 

我 们 系统 性 地 研究 了 IncRNA 在 蜜蜂 级 型 分 化 
表达， 找到 了 工蜂 和 蜂王 幼虫 发 育 早期 差异 表 
的 IncPRNA， 通 过 对 这 些 差异 表达 的 IncRNA 的 


P 


同 源 IncRNA 和 上 下 游 编码 基因 的 功能 分 析 ， 我 们 


为 ，IncRNA 的 确 参 与 了 蜜蜂 级 型 分 化 的 过 程 ， 
可 能 通过 影响 上 下 游 基 因 的 转录 和 功能 
行 方式 在 蜜蜂 幼虫 的 多 细胞 组 织 发 育 、 神 经 系统 
育 、 转 录 调 控 的 过 程 中 起 到 了 重要 的 调控 作用 . 
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The Function of IncRNAs in The Caste Determination of The Honeybee’ 
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Abstract Caste determination in the honeybee has been demonstrated to be decided by the royalactin in the royal 
jelly. And the differentially expressed coding genes in caste determination of the worker and queen larvae have 
also been well studied. Our experiments identified that IncRNAs are also differentially expressed in the caste 
determination of the larvae, which indicates that IncRNAs are involved in the caste determination. Detailed 
analysis showed that in the early development of the larvae, IncRNAs play a regulatory role in multicellular 
organismal development, nervous system development and regulation of transcription by affecting the expression 
and function of the nearby coding genes. 


Key words honeybee, caste determination, IncRNA 
DOI: 10.16476/j.pibb.2015.0110 


*This work was supported by grants from Chinese Academy of Science Strategic Project of Leading Science and Technology (XDA01020402) and 
Hi-Tech Research and Development Program of China (201244020402, 2012AA02A202). 

**These authors contributed equally to this work. 

*** Corresponding author. 

CHEN Sheng-Lu. E-mail: slchen@zju.edu.cn 

ZHANG Shao-Wu. E-mail: shaowu.zhang@anu.edu.au 

CHEN Run-Sheng. Tel: 86-10-64888543, E-mail: crs@sun5.ibp.ac.cn 

Received: April 22, 2015 Accepted: June 1, 2015 


